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ABSTRACT 

Neutron stars are the densest objects known in the Universe. Being the final 
product of stellar evolution, their internal composition and structure is rather 
poorly constrained by measurements. 

It is the purpose of this paper to put some constrains on the mass and mo- 
ment of inertia of neutron stars based on the interpretation of kHz quasi-periodic 
oscillations observed in low mass X-ray binaries. 

We use observations of high-frequency quasi-periodic observations (HF- 
QPOs) in low mass X-ray binaries (LMXBs) to look for the average mass and 
moment of inertia of neutron stars. This is done by applying our parametric 
resonance model to discriminate between slow and fast rotators. 

We fit our model to data from ten LMXBs for which HF-QPOs have been seen 
and the spin of the enclosed accreting neutron star is known. For a simplified 
analysis we assume that all neutron stars possess the same properties (same mass 
M* and same moment of inertia /*). We find an average mass M* p^ 2.0 — 2.2 M Q . 
The corresponding average moment of inertia is then i* ps 1 — 3 x 10 38 kgm 2 
0.5 — 1.5 (10 km) 2 M Q which equals to dimensionless spin parameter a ps 0.05 — 
0.15 for slow rotators (neutron stars with a spin frequency roughly about 300 Hz) 
respectively a ps 0.1 — 0.3 for fast rotators (neutron stars with the spin frequency 
roughly about 600 Hz). 

Subject headings: Accretion, accretion disks - Stars: neutron - Equation of state - 
Dense matter - Relativity - X-rays: binaries 
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Introduction 



Neu tron stars are excelle nt astrophysical laboratories to test matter above nuclear 



density ( jPage k, Reddy 



20061 ). Unfortunately, there is nowadays no way for nuclear 
physicists to investigate matter at such extremely high densities in laboratories. Moreover, 
because of the lack of knowledge about the behavior of particles in these extreme regimes, 
there is yet no consensus on a satisfactory equation of state for nucleons. Many modern 



equations of state have been propos ed, based on non-re 



help on relativistic field theory (see 



Haensel et al 



ativistic approximations or with 



(120071 ) and references therein). These 



equations of state at or above nuclear density predict different mass to radius relations 



for neutron stars. The answer or a piece of it cou 



laboratories but from the sky (ILattimer fc Prakash 



d maybe come not fro m terrestrial 



2007 



OzeletaL 



20101 ). Indeed, it has 



been claimed that measuring the mass and the radius of neut ron stars wil 



the p roposed equations of state and to reject some of them (IMiller et al. 



help to constra in 



1998 



Lattimer 



20071 ). 



HF-QPOs observations in LMXBs is a unique tool to test gravity in the strong field 
regime and to learn about the behavior of particles at high densities. Further detailed 
observations and modelling of QPOs for individual objects will help getting more insight 
into the properties of individual accreting neutron stars. 

How can we then estimate their mass and radius? In binary neu tron stars showing up 



as pulsars, the task is relatively easy (IThorsett fc Chakrabarty 



19991 ). The very accurate 



clock furnished by the pulsar serves a s an efficien t instrument to deduce the orbital motion 



20061 ). Such techniques have been successfully 



and other parameters in this system (INicd 
applied by numerous author s , findi ng masses aggregating around 1.4 M (for a summary, 



see e.g. 



Lattimer k Prakashl (120071 )1 



For neutron stars in LMXBs, the situation is less favorable although some attempts 
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have been made. For a rece nt review on d i fferent metho ds to comput e neut ron star 



parameters, see for instance 



Bhattacharyyal (120101 ) or also 



Zhang et al 



(120071 ). During 



their life, neutron stars in binaries accrete matter from their companion, an amount 
which can reach a substantial fraction of their initial mass. Not surprisingly, their final 



mass can deviate significantly from the fiduc ial 1.4 Mm. Actual 



stars seem plausible with M m 1.6 — 1.9 M (ICasares et al 



SAX J1808. 4-3658, such high masses were also found (iDeloye et al. 



y, high-mass neutron 



2006). For s ome pulsars like 



20081 ). 



Quasi-periodic oscillations (QPOs) seen in LMXBs can help to diagnose motion in 
strong gravitational fields and maybe solve the problem of determination of mass and radii. 
Although their estimates and related quantities are strongly mo del- dependent, picking up a 
particular QP O-model cert a inly h elps on mak i ng som e strong assertions about neutron star 



properties, see 



Miller et al. 



(119981 ) and 



Zhang! ( 120091 ) 



For a thorough review on X-ray variability and QPOs, see Ivan der Klisl (l2006bf ). 
Several models for high-frequency and low-frequency QPOs h ave been proposed as 



described in this review. Some invo ke resonance mechan isms (IKluzniak et al. 



other relativistic precession mo tion (IStella fc Vietri 



2007) 



2004 




Rezania & Samson 


2005) 


(Abramowicz et al. 


2007) 





19991 ) or MHD Alfven waves ( 



Zhang 



20051 ). However, some important problems are still unsolved 



In the present work, we show how to compute the average mass and moment of inertia 
of neutron stars by fitting kHz-QPO observations in LMXBs for slow and fast rotators. 
This paper is divided in two main parts. In Secf2] we brief ly sum marize the par ametric 



Petril (l2005af ) and in lPetril (l2005bf ). In SecEJ 



resonance model, more details can be found in 
we apply the model to ten LMXBs and deduce their key parameters. The conclusions are 
presented in SecJU Finally, Appendix |A] discusses the way to extend the model to allow for 
variable QPO frequencies, an important point with respect to observations. 
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2. Model and method 

In this section, we recall the main results of the model. The essential feature is the 
presence of a rotating neutron star which does not possess an axial symmetry about its 
rotation axis. The origin of the asymmetry can be due to the magnetic field tilted with 
respect to the rotation axis or due to anisotropic and inhomogeneous stellar interio r, 



Petri (12005ah has 



producing either a rotating asymmetric magnetic or gravitational field, 
shown that this induces some driven motion in the accretion disk due to a parametric 
resonance. Therefore the disk will show strong response to this drive by oscillating across 
the equatorial plane at some given preferred radii where the resonance condition is satisfied. 
More explicitly, remember that vertical resonance occurs whenever the vertical epicyclic 
frequency is equal to the perturbation frequency as measured in the locally corotating frame 

m\Sl{r,a*)-n*\ = 2 Kz ( r,a * ) (1) 

where m is the azimuthal number of the perturbation mode, fl(r, a*) the orbital frequency 
in the disk at radius r, the spin of the neutron star, a* a length related to the angular 
momentum J* = I* fi* by a* = J*/M*c, K z (r, a*) the vertical epicyclic frequency, J* the 
stellar moment of inertia and n an integer. The frequencies Q and k z are expressed for a 
test particle in Kerr space-time. They depend explicitly on the radius r and on the angular 
momentum a* as 



n(ra) = l — (2) 



d 0~ 

K z (r,a*) = ft(r,a,)^/l-4^ + 3^ (3) 

i?g = G M*/c 2 is the gravitational radius of the star, r — r/R g and a = a*/R g . We explicitly 
used the Kerr metric to find expressions for these orbital and epicyclic frequencies. However, 
this is not the best approximation for the exterior of realistic rotating neutron stars, since 
the quadrupole moment of the star usually causes large deviations from the gravitational 
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field that would be created by a sim pler Kerr b 



Kerr geometry was recently used by 



Torok et al 



a ck ho le geometry. Nevertheless, the 



(120101 ) to estimate the mass and the 



spin of neutron stars from the relativistic precession model. They argued that the metric 
well describes the exterior of rotating high-mass neutron stars and can be used when the 
non-rotating mass implied by the model is in the upper interval of masses allowed by 
the equations of state. This strongly supports our approach since already for a = the 
parametric resonance model implies a neutron star mass which is relatively high, around 
2M . 

Eq. ([3]) giving the vertica l epicy clic frequency in the Kerr approximation was first 



published 



Kato et al. 



Aliev fc Galtsov ( 



(1199a ) and 



1981). So me us eful properties are summarized in, e.g., 



Torok & Stuchlikl (120051 ) . 



From the known spin of the neutron star, we can deduce its angular moment by 

assuming a given value for the moment of inertia J*. Therefore, guessing a mass 
and a moment of inertia, we can solve quantitatively Eq. ([T]) for the orbital frequency Q 
and try to match observations of kHz-QPOs. 



For slowly rotating stars, o < 1, we retrieve the Newtonian expression 



tt(r, a*) « Q(r, 0) = « z (r, 0) = 
from which the solution of Eq. (OQ) follows immediately 



,3/2 



Q(r,0) 



mn 



mn±2 

The orbital frequency fl(r, 0) should remain smaller than this at the innermost stable 
circular orbit (ISCO) given in the non-rotating limit by 

/'u„, 2198 Hz I — ^ I = 1570 Hz 1 ' U ' 



(4) 



(5) 



(6) 



M,J \1AM Q 
This would give a first guess for the expected QPO frequencies, knowing the mass M*. 
Actually, because the spin frequency is well known from X-ray bursts for instance, we can 
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do better and include the angular momentum a into the description, but then the moment 
of inertia comes in as another free parameter. 

Several LMXBs have been observed with known spin rate and showing the twin peak 
QPO phenomenon. Depending on the neutron star rotation speed, they have been classified 
as slow rotator for v* < 400 Hz or as fast rotator for v„ > 400 Hz, {y* = Q*/2tt). For slow 
rotators, the twin kHz-QPO difference, Az/ obs w v*, is almost equal to the spin frequency 
while for fast rotators, it is equal to half of it, Ai> ohs v*/2. It is somet imes argued that 



this slow/fast rotator dichotomy is an artefact. iMendez fc Bellonil (120071 ) reexamined the 
data from all these sources and claimed that there is no clear trend in any segregation 
between them. They showed that the kHz QPO frequency difference Av is much more 
concentrated (mostly in the window [200,400] Hz, precisely a Gaussian with mean 308 Hz 
and standard deviation 36 Hz) than the range of neutron star spin (from 100 Hz to more 
than 600 Hz). Using a Kolmogorov-Smirnov statistical test, they found it highly improbable 
that Av an d are correlat ed. Av is almost constant and only weakly ^-dependent with a 



fit done by 



YinetaL 



((20071) who find 



< Av >f 



-0.19 u t + 389 Hz. 



(7) 



This makes the link betw een spin frequency and QPO frequency difference questionable. 



Mendez fc Bellonil (120071 ) went even further and made the strongest assumption of 
independence between both frequencies. Such hypothesis could rule out simple resonance 
models (those invoking linear oscillations and no inward motion of the flow for instance) as 
they claimed. However, this conclusion is not exactly true and would not hold anymore if 
some simplifying assumptions of any resonance model are left. Non-linear effects as well 
as a radially inward motion of the accretion disk can significantly change the oscillation 
frequency which becomes a function of the amplitude of oscillations and explicitly on 
time because the proper orbital, radial and vertical epicyclic frequencies vary when 
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matter approaches the neutron star surface. This was discarded so far. However, in this 
new extended picture, the spin frequency still plays an important role by triggering the 
resonance at some preferred radius, bringing the disk into off-plane oscillations that are 
slowly advected by the flow and drift downwards to the neutron star. Therefore, v* does 
not give a clear imprint to the precise kHz-QPO frequencies as it seems (not) seen in the 
data, but serves to launch the mechanism. Moreover, these motions occur due to matter 
flow influenced by gravity in a strong field regime, and thus the ISCO plays a central role. 
It is not the purpose of this paper to study the drifting and non-linear terms, which will 
deserve full attention in another work. Here, to give a taste, we only draw the basic lines of 
the consequences of these effects in Appendix [A] 

Any model predicting a fixed frequency ratio faces difficulties to explain the data 
since this frequency ratio is not only 3/2 or 4/3 (where m ost of the observations cluster), 



but covers a wider range as seen by 



Belloni et al 



( 120051 ) . Alt hough a strong line ar 



correlation exists, i t differ s significantly from the 3/2 ratio, see 



Abramowicz et al. 



Belloni et al. 



(2005b 



value first found by 



( 120051 ) and 



af). In addition, the frequency ratio clustering around the 3/2 



Abramowicz et al. 



(l2003al ) could be well explained by a uniform 
distribution of the lower and upper kHz QPO set in the source. The 3/2 peak in the 
observed ratio distribution comes from selection effects (sensitivity of measurement tools) 
since there is only a very narrow range of frequency ratio where bot h QPOs are sufficien tly 



stro ng in order to b e detec ted. The details can be found in works of 



Torok et al 



and 



Boutelier et al. 



fl2008aU bh 



(120101 ) who elaborated this issue. We emphasize that their results do 



not contradict the parametric resonance model. The question of the viability of su c h mod els 



Barret k Boutelier! (120081 ) 



remains fully open and subject to strong debates. Moreover, 
looked carefully at 4U1820-30 and found a gap of roughly 100 Hz in the QPO frequency 
distribution that is not attributed to selection effect and sharply peaked around a 4/3 ratio. 



For this special binary, it seems that some frequencies are disfavored. In other words, 
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within the orbital QPO interpretation, preferred radii exist within the disk, supporting the 
resonance model. 

When the star rotates slowly, its geometrized angular momentum a remains small and 
a first order expansion for fi SC o(5) with respect to a is possible. In the next section, we will 
show how to use this linear approximation to find severe constrains on the stellar mass and 
moment of inertia. 

Another way to tackle the resonance condition Eq. (pQ), in the general case for 
arbitrary a, is to work directly with the full expressions given by Eq. (T5])-(J3]). This requires 
a numerical algorithm to search for the allowed frequencies and is also done in the next 
section, Sec. [31 

A dozen LMXBs have been inventoried to exhibit the above mentioned behavior. The 
LMXBs sample used to fit our model for the kHz-QPO difference as measured by some 
other authors are summarized in Table [1] with appropriate references. We want our model 
to adjust to this set as close as possible by looking for appropriate mass and moment of 
inertia. Let us take an index % tracing this set of LMXBs by writing % e (LMXBs). For 
each binary, the observed twin kHz-QPO frequency difference is known as Az/° bs . Fixing M* 
and /*, we get a predicted Ai/ 4 model from our parametric resonance model. To evaluate the 
goodness of our fit, we introduce a merit function J 7 defined by summation over all the 
LMXBs and compare the discrepancy between predicted and measured QPO differences, 
such that 



with a statistical weight (jj. The summation should be understood over the set of observed 



the error in the observed QPO frequency difference for the binary labeled i. We use the 
Li-norm but other choices are possible like the L2-norm, although the latter being less 




8 



a obs/model 

systems. Aza 



are the observed/predicted HF-QPO frequency difference and cr^ 
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separation Av (in Hz) 


spin z/* (in Hz) 


ratio Avjv^ 


Millisecond pulsars 








XTE J1807-294 


179-247 


191 


0.94-1.29 


SAX J1808.4-3658 


195 


401 


0.49 


Atoll sources 








4U 1608-52 


224-327 


619 


0.36-0.53 


4U 1636-53 


217-329 


581 


0.37-0.57 


4U 1702-43 


333 


330 


1.01 


4U 1728-34 


271-359 


363 


0.75-0.99 


KS 1731-260 


266 


524 


0.51 


4U 1915-05 


290-353 


270 


1.07-1.31 


IGR J17191 


330 


294 


1.12 


SAX J1750.8-29 


317 


601 


0.53 



Table 1: The detail of the LMXBs with known twin kl 
source of the data, see for in stance 



Belloni et al. 



(120071 ) 



(120051 



z QPOs and spin frequenc ie s. For the 



20071 ): 



Zhang et al. 



(120061); 



Yin et al. 



van der Klis (12006a) and references therein. 
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robust. We found no significant changes when applying the second choice. In any case, the 
best fit corresponds to a minimum of the figure-of-merit function J 7 . Moreover, we tried 
other merit functions with no difference in the best fit parameters. 

Finally, some words about the mass-dependence on stellar rotation. Assuming the 
same mass as well as the same moment of inertia for all the set of neutron star binaries 
is a crude first guess. A detailed description of the inner structu re of rapidly rotating 



neutron stars is a diffi cult calculation only numerically treatable (jShapiro fc Teukolsky 



1983 



Camenzind 



20%. However, 



2007 ). For unifor m rotation, the mass increase is expected to be less than 



Lyford et al 



(120031 ) computed equilibrium configurations with differential 
rotation and found an increase up to 60%, even for moderate spin rate. The salient feature 
to keep in mind from all these studies is an increase of the gravitational mass with rotation. 

Thus, to better adjust the observations without handling all these complicated 
computations, we can however release the constant mass hypot hese and use a neutron star 



20071 ). The rotation 



spin dependent mass based on the following heuristic argument (IGhoshl ; 
of the neutron star, containing a fixed number of iV nucleons, increases its gravitational 
mass M(N, f2* 7^ 0) compared to the non rotating limit M(N, fi* = 0). Because kinetic 
energy is equivalent to mass and therefore induces gravitation, a simple relation between 
both gravitational masses is such that 



M(N, tt*)c 2 = M(N, 0)c 2 + l h Ql 



(9) 



For the remainder of the paper, we use lighter notations, setting M = M(N, 0) for the 
mass of a non-rotating neutron star and M* = M(N, fi*) for that of the same neutron star 
(i.e. equal number of baryons N) but rotating at an angular speed fl*. Eq. (jUJ) shows the 
quadrati c dependence on s pin the same functional dependence as the one from the 



study of 



Hartle fc Thornd ( 119681 ). The relative mass correction is therefore 



M M 2 <■- M. 



1.76 x 10 







1Q38 kg m 2 / V400 Hz 



(10) 
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So it remains small even for fast rotators. 



Results 



The above model and fitting technique is now applied to the dozen of fast and slow 
rotators. We emphasize that in order to make any prediction on the mass and moment 
of inertia, we have to take the same properties for the whole sample of accreting neutron 
stars in the observed LMXBs. Indeed, adopting different parameters for each system could 
significantly change the orbital and epicyclic frequencies. Most importantly, the frequency 
at the ISCO, scaling like 1/M* would change from one binary to another. But in our 
segregation between slow and fast rotators, the precise value of the orbital frequency at 
the ISCO is the salient feature to interpret the abrupt change in the twin peak frequency 
difference. A varying M* would shift this sharp transition to lower or higher frequencies 
from one binary to another. Thus the zero-th order choice, to highlight the general trend, 
is to keep the same mass for all neutron stars. 

Let us first give an estimate for the gravitational mass M* and moment of inertia I* 
along the following arguments. The geometrized spin parameter a is defined as 

-2 



0.145 



GMl V.10 38 kgmV V400 HzV \1AM Q 

From this expression it is clear that it remains small compared to unity, and this even 



11) 



for fast rotators. In this case, to first order in a , the orbital 



( iKluzniak fc Wagoner 



1985 



Kluzniak et al. 



1990 



Miller et al. 



requ ency at the ISCO is 



19981 ) 



z/ isco (5) = 2198 Hz (1 + 0.75 a) ^ 



(12) 



Assuming that a < 0.3, an inaccuracy introduced by this simplification with respect to 
the Kerr solution (due to neglecting higher order terms in j) is smaller than 10%. We put 



13 



explicitly the spin dependence through a on the left hand side for latter convenience. And 
therefore the relation between M* and I* becomes approximately 



2 GMl 



3 71 Cl>* 



- 1 



(13) 



2198 Hz J M Q 

According to our parametric resonance model, for slowly rotating stars, the twin kHz-QPOs 
are given by v\ = 2 and v\ = 3 v*, where the superscript s stands for slow. Because uf, v\ 
are interpreted as the frequencies of the orbital motion, they need to be less than that at 
the ISCO 



(14) 



For increasing spin of the neutron star u t , at some point, v\ will approach and eventually 
overtake v- xsco . Thus v\ will be forbidden as a HF-QPO. As a consequence, the next two 
dominant twin kHz-QPOs are identified as v\ = 1.5 and Therefore, the QPO 



frequency difference Av/v* = (v?. — v\ 



the data taken from 



Iv* jumps suddenly 



Mendez fc Bellonil ( 120071 ); Ivan der Klid ( 120081 ). this should happen in the 



rom 1.0 and 0.5. According to 



neutron star spin range G [363,401] Hz. This is probably the most salient feature in the 
slow against fast rotator discrepancies. Fitting these data requires that the switching from 
slow to fast rotator occurs for neutron star spin between 363 Hz and 401 Hz. More precisely, 
for u* < 363 Hz, Av/v* ~ 1 which we interpret as no effect on motion in the observable disk 
from the presence of an ISCO. This implies that z/ isco (363 Hz) > 3 = 1089 Hz, we put the 
spin rate into coma to distinguish between different rotators, an essential remark for our 
constrains. Next, for v* > 401 Hz, Avjv^ ~ 0.5 which we interpret as a clear signature of 
the ISCO. This implies that z/i SCO (401 Hz) < 3 v* = 1203 Hz. Express in terms of the ISCO, 
the transition from slow to fast rotator should happen when the two conditions below are 
satisfied 



z/ isco (363 Hz) > 1089 Hz 
z/ isco (401 Hz) < 1203 Hz. 



(15) 
(16) 
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This condition supplemented with the relation Eq. flT3|) sets two constrains on M* and /*, an 
allowed region in the (M*, /*) plane. Next, a third bound for the couple (M*, /*) is possible 
along the following lines. For fast rotators, the ISCO is clearly taken into account. But for 
the highest accreting system with = 619 Hz, the ratio is still Ai//z/* 0.5, the upper 
kHz-QPO being v\ = 2 = 1238 Hz and the lower kHz-QPO being v{ = 1.5 1/» = 929 Hz. 
We conclude that /or t/ws particular system 



The last and general constrain is that there is no naked singularity in the Kerr metric or 
stated mathematically, \a\ < 1. In terms of the moment of inertia, it means that 



The less favorable case (most restrictive one) corresponds to = 619 Hz. This leads to 



For later use, we introduce Jo = 10 38 kg m 2 . All these constrains, Eq. ( lT5|) -(|T8l). are 
summarized and shown in a (M*, JJ-plane depicted in Fig. [TJ The hashed regions are 
forbidden and only a small area in white survives around the first diagonal in the figure. 
This plot clearly emphasizes the existence of a lower and upper bound for both the mass and 
moment of inertia. We found the minimum values to be M* min = 1.9 M and i* m i n = 0.98 Iq 
whereas the maximum ones are M* max = 2.9 M Q and i* max = 19.1 Iq- Neutron star structure 
models predict 7* close to or slightly above I so that we will favor the lower bounds and 
expect masses in the vicinity of 1.9 M & . 

In a last step, we use the figure-of-merit function J 7 , Eq. (JH]), to fit the data and also 
the variation of mass with spin rate according to Eq. (Q. We span a vast range in the 
(M, .Z*)-plane to compute the merit function. Note however the subtil change in unknowns 



z/ isco (619 Hz) > 1238 Hz. 



(17) 




(19) 
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compared to the previous analytical study. Now we use a constant mass for the non-rotating 
limit, M, instead of a constant mass for the rotating star, M* = M + 7* f2 2 /2c 2 . Actually, 
the discrepancy between both approaches is small and can be neglected at the end of the 
study. 

The results presenting the isocontours of the merit function are summarized in Fig. [2J 
The resulting region for minimization of J 7 is marked in red and shapes very similarly to 
the previous diagram (M*, /*), Fig. [TJ The most probable mass and moment of inertia are 
M ss 2.0 - 2.2 M Q and h « 0.5 - 1.5 (10 km) 2 M . The best fit according to these values 
is shown in Fig. [3] where the spin rate is plotted on the x-axis and the twin kHz-QPOs 
difference normalized to the spin rate is plotted on the y-axis. First, we retrieve the 
segregation between slow and fast rotators at the correct frequency as expected. Next, for 
fast spinning stars, the theoretical curve agrees very well with observations. Nevertheless, 
for slow rotation rates, the spread around unity is significant and cannot be explained in 
a straightforward way by our predictions. Clearly, some refinement of the model is still 
needed and under investigations, including other aspects of the plasma flow around an 
accreting magnetized neutron star. 



4. Discussion and Conclusion 

In this paper, we investigated further the consequences of forced oscillations induced 
in an accretion disk to explain the twin kHz-QPOs in LMXBs. Our model is able to 



discriminate between slow and fast rotator as already shown in iPetril (j2005al ). Moreover 



with help on new data from a dozen rotators, we were able to constrain the average mass 
and moment of inertia of neutron stars. We found for the best fit M ~ 2.0 — 2.2 M & 
and 7* ~ 0.5 — 1.5 (10 km) 2 M & . Whereas the moment of inertia gives roughly the same 
value as those obtained from independent ways by solving the stellar structure with several 



- 16 - 



Constrain on M 4 and I* 




2 2.2 2.4 2.6 2 

M, /M 



Fig. 1. — The four constrains Eq. (|T5l)-(|T8l). in the (M#, J*)-plane, labeled with the spin 
frequency or marked by singularity for Eq. (fl8"|) . The hashed regions are forbidden. The 
minimum allowed mass and moment of inertia are around M ~ 1.93 M & and J* ~ 0.98 Jo- 



Merit function 




1.6 1 . 



2 

M/M 



2.2 2.4 



Fig. 2. — Isocontours of the merit function J 7 vs the fiducial mass of the non- rotating neutron 
star M, normalized to the solar mass M , and vs its moment of inertia /*, normalized to 
Iq = 10 38 kgm 2 . The minimum value of / lies around M ~ 2.0 — 2.2 M & and ~ 1 — 3/q- 
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equations of state ( iWorley et al. 



20081 ). the neutron star mass appears rather large. This 
effect could be an artefact of its constancy from one binary system to another. Better fits 
suggests to look at each system individually and remove the constant mass approximation 
for the whole set of LMXBs, leading to a spread in the mass distribution function for 
neutron stars. But this requires a much more detailed separate analysis of each binary 
with their own specificities (accretion rate, magnetic field strength for instance) and better 
observations. New time analyzing instruments like the HTRS (High Time Resolution 
Spectrometer) project on board IXO will give more insights into supra-nuclear matter and 



strong gravity physics (IBarret et al 



20081 ) 
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Slow vs fast rotator 



1 }} 






1 — ^ 







200 300 400 500 600 



Spin frequency v, (Hz) 

Fig. 3. — Observations of slow and fast rotators (red triangles) and fit obtained by our model 
(blue solid line). The best fit parameters are M « 2.0 — 2.2 M & and /* ~ 1 — 3 Iq. 
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Toward a more realistic model 



So far our linear parametric resonance model predicts fixed radii where the resonance 
conditions Eq. (JT|) are satisfied. Therefore the orbital frequencies remain also constant, 
leading to fixed kHz-QPOs. This can only be an approximation since oscillations are 
non-linear and gas or particles drift slowly towards the neutron star due to accretion. In 
other words, advection increases the orbital and vertical epicyclic frequencies and puts 
the system (particles) out of resonance. Actually, if non-linear terms are retained, the 
proper frequency of the vertical excursions depends on the amplitude of these oscillations. 
Therefore, it is possible that the excitation and proper frequencies adjust themself in such 
a way to maintain the oscillator in high amplitude motion. We call this a parametric 
auto-resonance mechanism. Its explanation is given in more details along the following 
lines. 

The idea of non-linear resonance to expl ain QPO observ a tions has already b een 



discussed by many autho rs, see for instance 



Abramowicz et al 



Rebuscol (120041 ): 



Horakl (12004 . 



20051 ): 



( I2003bl ). Although they considered a resonance between oscillatory 
modes different from those relevant to the model presented within this paper, their results 
show also how non-linear phenomena can drastically improve their models. 

Let us see how parametric auto-resonance works. Add a non-linear cubic term in the 
usual Mathieu equation governing the vertical displacement z(t). Note that, to first order, 
a quadratic term would not lead to a change in proper frequency with amplitude so a cubic 
term is more relevant for our discussion. Thus, our non-linear oscillator takes the form 
below 



d 2 z 

^2 + nl(t) [1 + h cos(m(fi(t) - ft*) t)]z = p z 3 



(Al) 



h is the strength of the excitation. Now, an important new fact is that the vertical 
epicyclic K z (t) and orbital Q(t) frequencies are time dependent. These variable coefficients 



-24- 



mimic the particle radial drift. How does the position of the resonant particles evolve with 
time due to loss of angular momentum expected from accretion? A simple argument to get 
the temporal dependence is the following. Assume that the thin accretion disk possesses a 
power-law axisymmetric surface density S(r) such that 

£(r) = £ (£f (A2) 

where a is the power law index and E the surface density at r . The flow starts to accrete 
at an initial speed v (directed radially inwards such that v = —v e r ) at radius r . The 
conservation of mass implies 

2 7rrS(r)i;(r) = 2 7rr S Wo (A3) 

Integration with respect to time, using the fact that for a test particle v(r) = —dr/dt 
(projection along — e r ), we get for 

r(f) = [r a+2 - (a + 2) r " +1 v (t - *„)] ^ (A4) 

where tq = r(to) is the initial position of the particle. Specializing to a uniform density disk 
model, i.e. a = 0, the particle falls onto the neutron star along the trajectory 

r(t) = ^Jr 2 -2r v (t - t ) (A5) 

For convenience, in the remainder of this paper, we will use this expression for the particle 
radial path. Furthermore, we make a shift of time by the replacement (t — to) — > t. This 
temporal dependence on radius r(t) governs the time evolution of the orbital Q(t) as well 
as the vertical epicyclic K z (t) frequencies. Parametric auto-resonance will settle in if the 
excitation is efficient enough to maintain phase locking. In so me special cases, we can 



look for semi-analytical solutions by the method of averaging ( iKhain fc Meerson 



200l|). 



The underlying idea is to smooth the fastest time scale Q(t) and to keep track of only the 
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secular amplitude change and phase evolution of the oscillation, compared to an harmonic 
oscillator. 

We apply this technique to our resonance model. First, for small enough times 
v t <C r , an expansion of radius and frequencies yields (assuming a non-rotating body, 
5 = 0) 



r(t) ~ r — v 1 



1 3 f o t 
2 r 



fi 

K z (t) 



IGM 



(A6) 
(A7) 

(A8) 
(A9) 



The non-linear parametric resonance model for vertical motion becomes 

2 



dP + V 2 r 



1 + h cos { m ( fi ( 1 + 7^ _ ~ — 



t 



z = /3z 3 



(A10) 



Next, the method of averaging looks for solutions expanded into z(t) = A(t) cosi9(t), the 
amplitude being A(t) and the phase being $(t), A(t) varying on a timescale much longer 
than the period of oscillations. It is preferable to introduce a new phase defined by 



if) = - fl 1 



(All) 



to get rid of the fastest timescale represented by Qq. Let us have a look on the behavior of 
the first resonance. Specializing to this particular case corresponding to n = 1 and m = 1, 
the resonance condition, in the Newtonian limit, is Qq = fi*/3. Inserting into Eq. 
the vertical motion satisfies 



d 2 z Ql / 3v t 

1 — *(H — 

dt 2 9 V 2 r 



(A12) 



The straightforward way is to solve numerically this second order differential equation with 
appropriate initial conditions. We show an example of numerical integration of Eq. flA12j) 
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in Fig. HJwhen the particle enter in the first resonance, n = m = 1. The full solution, z(t), 
is not plotted because of the small timescale, whereas the evolution of the amplitude is 
shown by a solid blue line, |z(i)|. A clear increase of the amplitude with respect to time is 
demonstrated. This shows that resonance can occur even for variable orbital and excitation 
frequencies. Averaging provides another mean to solve approximately for the amplitude A 
and the phase if}. For accretion timescales much longer than the orbital motion, which is 
usually the case, these new unknowns evolve according to 

dA hAQ* sin(2^) 

~dt ~ 12 ^ ' 

dip htt*cos(2ip) 9/3 A 2 ft* u * 

H ~ 12 8aT + T^T 

A stationary solution corresponds mathematically speaking to a fixed point satisfying 

= ^ = 0. Inverting this system for the phase and amplitude, we find the real positive 
amplitude by 

sin(2^) = (A15) 

This analytical solution, solid red line in Fig. HI agrees very well with the straightforward 
numerical integration of Eq. (1A12|) . From that, we conclude that oscillations grow slowly 
as time goes. Simultaneously, the orbital period increases according to Eq. ( 1A7|) and 
consequently so does the kHz-QPO frequency. We recall that this describes the short time 
evolution of the particle just after entering resonance. In this limit, the fractional increase 
in the frequency remains small, a few percent. Other resonance, with different m,n will 
follow the same trend. Accretion and non-linearities allow the QPO peaks to drift towards 
higher frequencies and release the fixed frequency ratio result. 

The analysis of its long term evolution (with the full non-linear time dependence 
retained) requires more investigation and is left for future work. This discussion aimed 
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at emphasizing that the parametric resonance model extended to more realistic situations 
met in accretion disks (non linearity and chirped frequency) can explain an increase in the 
observed kHz-QPO frequencies by parametric auto-resonance. 



To conclude, representative refer ences about t h e real b eh aviour o 



and coherence times can be found in 
refer to 



T5rokl (120091 ) and 



] 3 arret et al. 



Horak et al 



(12005a 



Mendez 



;he Q PO amplitudes 



(120061 ) . Finally, we 



(120091 ) for a discussion on the relation between the 



strength of the twin QPOs and its possible explanation within the framework of a non-linear 
resonance theory. 

The study done in this appendix is very preliminary and will be included in an 
extended model of our parametric resonance. 
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h = 0.01, il = 10~ 6 




Fig. 4. — An example of growing oscillation amplitude solid blue line, due to the 

parametric auto-resonance. The parameters are h = 10 -2 , (3 = 10~ 2 and \l = Vq/tq — 10 -6 . 
The orbital frequency is normalized to unity 17 = 1 or equivalently = 3. The short orbital 
period of 2 7r is not plotted. The solid red line shows the quasi-stationary fixed point solution 
for the amplitude A(t), Eq. ( 1A16I) . in good agreement with the true solution. 



